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1. Introduction 

There is a longstanding discrepancy between experimental results for the like-sign dimuon 
asymmetry measured by the DO Collaboration [ 1 ]- [4] and the corresponding standard model pre¬ 
dictions [5]-[9] . In [10] an interesting connection between the dimuon asymmetry and the decay 
rate difference of neutral B d mesons was suggested: the measured enhancement of the dimuon 
asymmetry could also be explained by an enhanced decay rate difference of the neutral B d mesons. 
Moreover, A12/, is currently only weakly constrained by direct measurements. This was the mo¬ 
tivation for the study in [11] where indirect experimental constraints on possible new physics en¬ 
hancements of A T c i were studied; we present here the main findings of [11]. 


2. Neutral B mixing 

Due to electroweak interactions the neutral states B d and B d oscillate into each other, the time 
evolution of this system is given by solving the Schrodinger-like equation 


■ d (\B d )\ d (\B d )\ 
dt ^ ) 


where 


L d 



( M uMil 

\M d * - Mfj 



(2.1) 


Diagonalizing L d gives the physical eigenstates {\Bu >,|Bi >} with the masses M‘h.M d and the 
decay rates F^.P/. To provide a mathematical description of the mixing process it is useful to 
define the observables A M d = M'h — M‘[ and AT^ = T d H — V d . Theoretically A M d and AF^ can be 
calculated from the components of L d according to the formulas 


NM C , k, 21 m( 2 | and AT d ~ 2 j rf 2 1 cos ( (j ) d ) where <f> d = arg( - —p-). 


M c 


yd 

1 12 
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Only A M d and Ar^ are directly accessible in experiment, the phase (j> d can be calculated from the 
measurement of the semileptonic asymmetry 


4 = 


r d 

1 12 


Mf :2 


sin((f) d ). 


(2.3) 


By combining data from Belle, BABAR, DO, DELPHI and LHCb the following direct experimental 
bound for AT d is available 


A T” FAG 


T d 


( 0.1 ± 1 . 0 )% [ 12 ] 


which can be compared to the standard model prediction 


(2.4) 


AT sm 

— P- = (0.42 ±0.08)%[5]. 
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3. Sources of enhancement for AT d 

3.1 CKM unitarity violations 

Let X'l = V* d V q b, from the unitarity of the 3x3 CKM matrix the following condition is satisfied 
X'j + X‘' + X? = 0. However in different extensions of the standard model the previous equality is 
broken according to X£ + A/ + X’/ = 8f :KM \ for instance in 4th family studies 8q KM ~ 10 2 leading 
to a potential enhancement factor on AT,/ of 300%. 

3.2 Current-current standard model operators 

In this section we consider the possibility of having new physics effects in tree level standard 
model operators and its consistency with recent experimental observations, a similar study and its 
implications over the CKM phase y has been discussed in [13]. 

We are concerned with the following effective Hamiltonian 


^current . I AS I = 1 


4 Gp 


I v; d v p , b £ cf p '^)Qf p, +h.c. 

n 1 — ii r i=1 9 


(3.1) 


with 

Q p / = (d a pP) v . A (p' p b a ) v . A , 

Q P / = ( d a p a )v-A(p' P bP) V - A • (3.2) 

New physics effects will produce shifts over the standard model Wilson coefficients according to 


c i P 2 = c ia +ACf^. (3.3) 

To constrain the values of ACf p 2 we consider different observables depending on the up type quark 
structure indicated by the labels pp' in Eqns. (3.1) and (3.2). The most important bounds arise from 
the operators related with the transitions b —> uud, b —y cud and b —> ccd with Wilson coefficients 
C uu , C uc and C cc respectively. The decays B —> nn, B —? up and B —> pp provide limits for AC'"'. 
To calculate the bounds on AC" C the relevant processes are B —> D* + n and B° —r D^°h° with 
h° = 7t°. rj or co. Regarding that the decay width of the B d meson is dominated by the transition 
b -o- cud an extra constraint for AC'"' can be imposed from the decay rate itself. The shift AC CC 
is bounded by the process B X d y. the observable sin{ 2/3) and the semileptonic asymmetry a d 
defined in Eqn. (2.3). We visualize graphically the parameter space available for new physics by 
plotting the real and imaginary parts of ACj^, as an example we present the results for AC‘ 2 ,: in the 
left panel of Fig. (1). 

Finally the effects of Aover ATj can be estimated by using Eqns. (2.2) and (3.3) together with 
the expression 
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Re ACf(M w ) 



-4-20 2 

Re ACf(M w ) 


Figure 1: Left panel: Allowed parameter space in the Re AC^-ImAC^ plane. The constraints shown 
correspond to B —> Xjj (Blue), a d sl (Green) and sin( 2/3) (Red). The region contained within the dashed black 
line represent the combined constraints from the different observables. Right panel: contours of AT^/AT™ 
along with the combined constraints. 
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(0) |B rf > 


= - (Cf,cf) +2A c A ; ,r"^(Cf,Cf) + A 2 r;"’ rf (cr,Cf)) . (3.4) 


The enhancement factors can be read by plotting the contour lines for AT f / /AT s ( j^, we provide 
an example of this strategy in the right panel of Fig. (1). Our calculations show that there is 
room for new physics in C™ and C“ c , leading to enhancement factors of 1.6 in both cases. The 
most interesting result arises from CF' where an enhancement factor of 14 is allowed by current 
experimental results. All the bounds were calculated at 90% C.L.. 

3.3 Operators (db)(zT) 

In the third part of our analysis we consider the effective operators (db)( tt), they are well 
suppressed within the standard model, however they are not quite constrained by the experimental 
data available nowadays; the approach followed in this section is analogous to the study performed 
for B s mesons in [14]. Flere we take into account the following effective Hamiltonian 




4 Gf 


ij 


and study the scalar, vector and tensor Dirac structures of {db){ tt) 


(3.5) 
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Qsab = ( dP A b ) (zP B z) , 

Qv,ab = (df-P A b) (zy^P B z) , (3.6) 

Qta = (da^ v P A b) (za^ v P A z) 

with Wilson coefficients Csab, Cv A b and Cj.a respectively. Following the notation of [14]: A.B = 
/?,LandiV = (l±7 5 )/2. 

To get bounds on the Wilson Coefficients we use direct and indirect constraints. Direct bounds 
arise from particle transitions where the chain b —> dz + z~ appears at tree level; indirect estima¬ 
tions are derived from processes where the terms given in Eqns. (3.6) arise as a consequence 
of operator mixing, loop level corrections or both. The direct category includes the processes 
B c [ — y t + t~, B — > X d z 1 z and B + -A- n z l z . In the first case the following experimental bound 
is available Br(B d -A T ! T ) < 4.1 x 10 3 [15]. In the second and third cases currently there are 
not experimental measurements for the branching ratios. However we can estimate an upper limit 
for these observables by considering the space available for new B d decays, i.e. those that have not 
been measured yet. To fulfill this task we compare the ratio ^ — I j against ^ — I j and 

ignore new physics effects in the B s sector; our result is Br(B d -AX new ) < 1.5%. From B ( \ —y z z 
we get chirality independent bounds for Q, and Cy, whereas from B -A X d z + z~ and B -a n"z z 
we get constraints for all the Dirac structures shown in Eqns. (3.6). Indirect bounds are calculated 
from the branching ratios of the processes B,j -a X d y and B -a 7 T /./ /J . In the case of B d -a X d y 

the tensor element in Eqns. (3.6) mixes with the operator mediating the transition b -A dy. For 

the process B + -A 7t p ii there are contributions from mixing between the vector and tensor 
components of Eqns. (3.6) and the operators responsible for the chains b -A dy and b -A dl + l~. 

To quantify the effects of our ( db){zz) operators over Ar^ we first parametrize the new physics 
contributions in Fp through the factor A d 

r ?2 = rff%. (3.7) 

The relationship between A d and the Wilson coefficients of the operators in Eqns. (3.6) is estab¬ 
lished by the following inequalities 

\Ad\sAB < 1 + (0.4llQ;Qg)|Cs,As(/7h))| 2 

\Ad\vAB < 1 + (0.421q;q8) \CvAB(mb)\ 2 (3.8) 

\Aj\ta < l + (0.42+^ 3 g )|C^K)| 2 

here we are assuming single operator dominance; i.e. we are analyzing the effect of only one op¬ 
erator from Eqns. (3.6) at a time by switching on the corresponding Wilson coefficient and setting 
the others to zero. The subindex of |Aj| in Eqns. (3.8) indicates the operator under study. 

We finally arrive to the following enhancement factors at 90% C.L. 
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Figure 2: Bounds on the enhancement factor for AT</ as a function of the improvement on the experimental 
limits over the branching ratios for the processes Bj —> T~Z (purple), Bj —» X r /T 1 T ~ (blue) and Bj —> 
7T 1 T 1 T (red). Only the effect induced by the vectorial version of the (db){ft) operators is shown at 90% 
C. L.. 


< 1-6 

\&d\v,AB < 3.7 (3.9) 

< 1 - 2 . 

Stronger bounds for |A^| can be established if upper limits for the branching ratios of the processes 
B c i H> T 1 T , B — > XiiT + z~ and B + — > 71 1 T ‘ T are reduced experimentally; this effect is shown 
explicitly in Fig. (2) for the vector structure of the {db){ ft) operators. 

4. Like sign dimuon asymmetry 

Experimentally the DO collaboration has measured the raw like sign dimuon charge asymmetry 
A = Acp + Ahkg [l]-[4], after removing the backgrounds A^ g the remaining component Acp is the 
result of CP violation arising from neutral B meson interactions. Before 2013 it was assumed that 
Acp was caused only by CP violation in mixing between the B c j and B t j states leading to 


Acp x A d sl = C d a d sl + C s a s sl . (4.1) 

d s 

The terms a j are the semileptonic asymmetries and C,i s are proportionality constants depending 
on the production fractions, oscillation parameters and decay widths of B c j )S - Based on this inter¬ 
pretation DO reported in 2011 a measurement for A d sl that was in disagreement with the standard 
model prediction by 3.9a [3]. 
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Recently the contribution of CP violation in interference between decays with and without mix¬ 
ing was included in the analysis of App. According to [10] this implies the replacement of Eqn. 
(4.1) by 



(4.2) 


AL, 

r d ■ 


here 4D j s highly suppressed by the constant C'\ t in comparison with the contribution due to 


In principle Eqn. (4.2) may suggest that the tension between theory and experiment for App could 
be eliminated if Ar^ gets enhanced through non standard model physics. However a more detailed 
study reveals that Eqn. (4.2) is the first approximation towards the inclusion of CP violation in 
interference; instead of having a complete dependence on AT j it is expected that the relevant con¬ 
tribution to the interference depends on the components Tp and Tp, see Eqn. (3.4), in a more 
convoluted way than the one given in AFj [16]; this possibility is currently under further investiga¬ 
tion [17], 

5. Conclusions 

We have explored the possibility of having new physics effects on Ar t / within the Heavy Quark 
Expansion. Firstly we found that by breaking the unitarity of the CKM matrix by 10 2 we can 
have a deviation of 300% with respect to the standard model expectation on AF,/. Our analysis of 
dimension six tree level standard model effective operators shows interesting deviations over Ar^ M , 
the most remarkable example is given by ( dc)(cb ) where an enhancement factor of 14 is allowed 
by current experimental data. Finally in the case of the effective operators (dfi)(fr) we have found 
that AT C , can be nearly 4 times bigger than in the standard model. We mentioned that as a first 
approximation there is a relationship between Ar,/ and the dimuon asymmetry App measured by 
DO, describing briefly possible corrections that could be useful in reducing the gap between theory 
and experiment for App. Thus we strongly motivate some further experimental studies of Ar,/. 
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